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Abstract 



We study the two-flavor quark-meson (QM) model with the functional renormalization group (FRG) to describe the effects of 
. . . collective mesonic fluctuations on the phase diagram of QCD at finite baryon and isospin chemical potentials, hb and /i/. With only 
isospin chemical potential there is a precise equivalence between the competing dynamics of chiral versus pion condensation and 
LiJ that of collective mesonic and baryonic fluctuations in the quark-meson-diquark model for two-color QCD at finite baryon chemical 
potential. Here, finite - 3/U introduces an additional dimension to the phase diagram as compared to two-color QCD, however. 
'Ai zero temperature, the (yU/, //)-plane of this phase diagram is strongly constrained by the "Silver Blaze problem." In particular, the 
^ onset of pion condensation must occur at /// — m„/2, independent of fi as long as jj. + jj.; stays below the constituent quark mass of 
'the QM model or the liquid-gas transition line of nuclear matter in QCD. In order to maintain this relation beyond mean field it is 
crucial to compute the pion mass from its timelike correlator with the FRG in a consistent way. 
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Q^l. Introduction 

(D. 

• The phase diagram of Quantum Chromodynamics (QCD) 
.continues to receive enormous attention both experimentally 
and theoretically worldwide [1-3]. QCD is expected to re- 
veal a rich phase structure in the plane of temperature T and 
baryon chemical potential fiB- Monte-Carlo simulations in lat- 
tice gauge theory provide a powerful non-perturbative first prin- 
ciple approach to QCD. At finite baryon chemical potential, 
however, they are restricted by the fermion-sign problem, the 
f-*) "fermion determinant becomes complex, in general, and it can 
p\| .therefore no longer be interpreted as part of a probability mea- 
"sure. There are various ways to deal with the sign problem [4], 
^ .but the region of large baryon chemical potential and low tem- 
"perature, fiB '» T, essentially remains inaccessible to lattice 
rN .simulations. This is one motivation to study effective models 
with QCD symmetries such as (Polyakov-)Nambu-Jona-Lasinio 
((P)NJL) [5-7] or (Polyakov-)quark-meson((P)QM) models [8- 
11], in order to describe the expected gross features of the QCD 
phase diagram at finite baryon chemical potential. 

The sign problem also motivates studies of QCD-like the- 
ories without this problem such as QCD with isospin chem- 
ical potential yU/ or two-color QCD [12, 13]. In both cases 
the fermion determinant remains real, and for an even num- 
ber of degenerate quark flavors positive at finite density. Con- 
sequently, lattice studies can be performed for both, QCD at 
finite isospin density [14—17] as well as two-color QCD at fi- 
nite baryon density [18-21]. In both these cases one studies 
a bosonic superfluid with BEC-BCS crossover towards higher 
densities in the low temperature, finite density region, however. 
An interesting QCD replacement with fermionic baryons is the 
exceptional G2 gauge theory [22] which can also be simulated 
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at finite density without sign problem [23] . Such lattice simu- 
lations allow to benchmark and cross-check the effective model 
calculations [24] as well as non-perturbative functional contin- 
uum methods for full QCD in general [25, 26]. 

In this letter we use the functional renormalization group 
(FRG) to study collective mesonic fluctuations in the two-flavor 
QM model for QCD at finite baryon and isospin chemical po- 
tential. As compared to the standard (fiB, T) phase diagram, 
finite isospin chemical potential thereby induces an imbalance 
between up and down quarks as in neutron stars or heavy ion 
collisions, for example. The region of high baryon density, 
where this becomes relevant, cannot be described by a QM 
model without any baryonic degrees of freedom, of course. On 
the other hand, in color superconducting quark matter at even 
higher density, isospin chemical potential will eventually de- 
stroy BCS pairing probably leading through an FFLO phase 
into an unpaired state [27, 28]. The other way round, finite 
baryon or isosymmetric quark chemical potential ji- fiBl'i will 
lead to an imbalance of the Fermi spheres of up and anti-down 
quarks and thus destroy pion condensation in the (yU/, T) phase 
diagram, transgressing through an FFLO phase, likewise [12]. 

Moreover, the standard QCD phase diagram with baryon 
chemical potential is related to that with isospin chemical po- 
tential outside the pion condensation phase via orbifold equiva- 
lence in the large A^^ -limit [29]. Just as two-color QCD [13], the 
latter can be described at sufficiently low temperatures within 
the framework of chiral effective field theory [12, 30, 31] and 
random matrix theory [32, 33]. At zero temperature, chiral per- 
turbation theory predicts a second order quantum phase transi- 
tion at yuj = m„l2, the onset of pion condensation. This is an 
exact property whose explicit verification from the grand parti- 
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tion function is called the Silver Blaze problem [34]. It relates 
the onset of pion condensation to the (vacuum) pion mass m„. 
The charged pions couple to 2yU/, so their excitation thresholds 
are given by m„ ± 2/i/, and the zero-temperature limit of the par- 
tition function must remain independent of fXi below the critical 
value of the isospin chemical potential at which m„ — 2/j.j = 0. 
The two-flavor theory at /// = is invariant under the full 0(3) 
flavor rotations in no,n+,7T-. Isospin chemical potential is as- 
sociated with the conserved charge corresponding to the 0(2) 
subgroup of rotations in the n+,n- plane. It acts as an exter- 
nal field which explicitly breaks the 0(3) down to 0(2), but 
the zero temperature limit of the partition function remains un- 
changed as long as its strength is not suflicient to excite charged 
pions. When the isospin chemical potential reaches its criti- 
cal value at /ij - m„/2, on the other hand, the resulting pionic 
Nambu-Goldstone excitations of zero energy lead to the sponta- 
neous breaking of the corresponding 0(2) symmetry, and Bose- 
Einstein condensation (BEC) of charged pions occurs. 

Considering non-zero values for both isospin and isosym- 
metric quark chemical potential, more general constraints arise 
from the Silver Blaze property in the zero temperature (;U/,/i) 
phase diagram as we will discuss below. We furthermore present 
the phase diagram of the QM model for two-flavor QCD in the 
three-dimensional parameter space of temperature, isospin and 
isosymmetric quark chemical potential. As compared to previ- 
ous model studies, mainly within the NJL model [35^2], the 
ERG approach is thereby well equipped to include the corre- 
sponding collective mesonic fluctuations in the descripiton of 
chiral versus pion condensation. 



2. Quark-meson model with isospin chemical potential 

With the standard two-flavor chiral symmetry-breaking pat- 
tern, S U(2)i X S U(2)r SU(2)v, finite isospin chemical po- 
tential leads to a further explicit breaking of the S U(2)v isospin 
symmetry down to the rotations about the iso-three direction, 
i.e., on the quark level to SU(2)v U(l)y\ In the pion con- 
densation phase at ii[ > ii'j(T), with /ij(0) = m„/2, this rem- 
nant chiral symmetry then spontaneously breaks down to a dis- 
crete Z2 synnmetry, U(l)y^ Z2, with a non- vanishing pion- 
condensate as the order parameter and one Nambu-Goldstone 
boson. As the isospin chemical potential is further increased, 
the pion condensate increases and the chiral condensate de- 
creases corresponding to a rotation of the vacuum alignment 
in a BEC-BCS crossover. The approximate chiral symmetry 
gets partially restored, and the asymptotic Z2 x J/(l)^' symme- 
try for /i/ — > 00 in the pion condensation phase agrees with that 
obtained in a hypothetical world with isospin chemical poten- 
tial but without any chiral symmetry breaking Ot'SB) and hence 
quark mass = as summarized in Table 1 . 

In this letter we use the QM model as a chiral effective 
model in which this symmetry pattern is realized. The Eu- 
clidean Lagrangian of the QM model at finite quark chemical 
potential = but zero isospin chenucal potential is given 
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Table 1: Symmetry breaking patterns in the charged pion condensation phase, 
by 

-Cqm = ijf(^ + g(cr + J'rs^^ - Wo) ^ 

+ \(dt^(rf + 1(5^70" + U(o-^ + 7f2) - ca, 

with mesonic potential U(cr^ + #) = ^(o"^ + t?^) + ^(o"^ + ^^)^ 
and explicit;k'SB term ccr corresponding to a finite current quark 
mass. The Lagrangian is invariant under 0(3)-isospin rotations. 
Eor the iso-three direction these are of the form 

ij/ e'^'^i^, il/' ^ e~"^^\lr\ n± = nx±in2 ^ e^'^'^Tii, (2) 

with an associated conserved current 

JI = i'k'^iyn^ + 2i(7T-dft7:+ - n+d^n-). (3) 

Isospin chemical potential is now introduced in the usual way, 
by adding a term niQ^ with the associated conserved charge to 
the corresponding Hamiltonian. Following the standard deriva- 
tion [43], one obtains the Lagrangian with finite isospin and 
isosymmetric quark chemical potential, conveniently written in 
the basis of up and down quarks with i/r = (^^r^, il/^)^ as 

XgM+w =>I'S'q4' + \idi,af + \(di,nof + U(p\d^) - ca 

I (4) 
+ - {(d^ + 2ixj^^)n^(d^ - 2/i,j0);r_) , 

where = cr^ + tTq (with ttq = tt"), cP- = tt+tt-, and 

Setting n = temporarily, it is evident in this representation 
that there is a precise map between the quark-meson model for 
QCD with isospin chemical potential ni and the correspond- 
ing quark-meson-diquark model for two-color QCD [44] with 
baryon chemical potential provided by the following modifica- 
tions: Instead of the bi-spinors in flavor here, one introduces 
bi-spinors in color consisting of say a red quark xj/r and a green 
charge-conjugated antiquark ij/'^ = tiCi/'J, with n for complex 
conjugation in the flavor 5 1/(2). The charged pions tt- and 7:+ 
are then replaced by a scalar diquark-antidiquarkpair, A and A*, 
coupled to the baryon chemical potential fis = 2/1 for A'c = 2 in- 
stead of 2/1/ for the charged pions in QCD with isospin chemical 
potential, and the neuttal pion ttq is replaced by the isovector of 
three pions ^, with no — » fi^ in 5 r ^ which increases the number 
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of would-be-Goldstone bosons for^SB from 3 to 5 in two-color 
QCD with two flavors. This identification is vahd only for the 
matter sector, of course, as the gauge sectors of both theories 
are fundamentally different. 

3. FRG flow equations 

3.1. Effective potential 

The functional renormalization group is a well-established 
non-perturbative tool to study critical phenomena in quantum 
field theory and statistical physics. It describes the evolution 
of the scale-dependent effective average action Fk from a mi- 
croscopic bare action specified at a UV cutoff scale Auv to the 
full quantum effective action for — > in terms of a functional 
differential equation [45] 
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Figure 1 : Diagrammatic representation of the flow equation for the effective ac- 
tion. Solid (dashed) lines represent full field and RG scale k dependent fermion 
(boson) propagators and the circles insertions of d,Rk(q)- 



Ansatz for the scale-dependent effective action is chosen to com- 
ply with the expected symmetry-breaking pattern. In particular, 
because isospin chemical potential explicitly breaks the (9(4) 
flavor symmetry of the mesonic sector down to (9(2) x (9(2), the 
effective potential Uk must be allowed to depend on two invari- 
ants, - cP- + tTq and cp- - 71^^ + ji^ - n+n^. Therefore, our 
Ansatz for the effective action is given by 



{n,m) 



denotes the nth functional derivative with respect Tk[p^,<f'] - T ^ I d^x -Cqm+i 



where 

to fermionic and m th functional derivative with respect to bo- 
sonic fields of the effective average action and the trace is taken 
over internal degrees of freedom and momentum space. Fig. 1 
shows a diagrammatic representation of Eq. (6). Although the 
flow equation has one loop structure, it is exact. Apart from 
the general requirements to act as the desked infrared regula- 
tor consistent with the boundary conditions [46, 47], there is a 
certain freedom in the precise choice of the fermionic/bosonic 
regulator functions Rkp/RkB which can be exploited to minimize 
truncation effects. Here we use 
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with X.QM+11, from Eq. (4). At the UV cutoff scale we assume 
an 0(4) symmetric potential of the form 
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Inserting the Ansatz (9) into the Wetterich equation (6), the flow 
equation for the effective potential can then be derived entirely 
analogously to that of the quark-meson-diquark (QMD) model 
for two-color QCD [44], yielding. 
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which are three-momentum analogues of the optimized regula- 
tor for the local potential approximation (LEA) [48] . They have 
the particular advantage that the Matsubara sums in the flow 
can be evaluated analytically. The regulators suppress the prop- 
agation of momentum modes below the scale k. Quantum fluc- 
tuations from all momentum modes are included by integrating 
the flow equation (6) from the bare classical action at Auv down 
to ^ = 0. From n functional derivatives of (6) one obtains the 
flow equations for the n-point vertex functions which in turn 
contain up to (n + 2)-point functions. Just as Dyson-Schwinger 
equations, the flow equations for all n-point functions form an 
infinite hierarchy which requires truncations to obtain a closed 
system of equations suitable for non-perturbative solutions. 

One frequently used truncation scheme is the derivative ex- 
pansion, which has been applied successfully to 0{N) and quark- 
meson models. In contrast to the NJL model, the QM model 
contains explicit mesonic degrees of freedom in the Lagrangian 
and is well suited to study mesonic fluctuations with the FRG. 
At leading order in the derivative expansion a scale-dependent 
effective potential Uk is obtained without wave function renor- 
malizations and with scale-independent Yukawa couplings. The 
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Here, E± denote A;-dependent quark excitation energies, 

£± = ^f^27(£~7^, with E, = ^k^+gW, (12) 
and for the neutral pion, we introduced 



E^ = V^272t7^, 



(13) 



with U,- = ^,Ud= 1^, and f/rr = etc. used as short- 

hand notations for the derivatives of the scale-dependent effec- 
tive potential here and below. The w, are the positive square- 
roots of the three poles in a? - -p^^ of the scale-dependent 
boson propagator in the 3x3 subspace of n^, and cr mesons 
which mix at finite isospin density. They are the roots of the 
same cubic polynomial as given for the QMD model in [24, 44], 
and the R, are the corresponding residues. For d^ = n+n^ - 
0, without pion condensation and pionic fluctuations, one has 
a>o - Eg- - ^jk^ + 2Ur + 4p2f/,.,., = 1 for the sigma meson, 
and wi 2 - Ejj + Ifj^s with R\ 2 - cl)\^2/E7t for the charged pions. 
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For zero isosymmetric quark or baryon chemical potential, 
fi - 0, Eq. (11) is identical to the corresponding flow equa- 
tion for the eff'ective potential of the QMD model for two-color 
QCD at finite baryon density [44] with Nc — 2, fxi ^ fi, and 
a multiplicity of 3 for the first term on the right in Eq. (11) to 
account for the usual three pions replacing the neutral one here. 

3.2. Pion pole mass and flow of the 2-point function 

It is common practice in quark-meson model studies to fix 
the parameters by adjusting the screening masses, i.e., the eigen- 
values of the Hessian of the eff'ective potential at its minimum, 
to the physical ones such as the pion and sigma meson masses 
here. With Bose-Einstein condensation of diquarks in two-color 
QCD or charged pions here, however, the zero temperature on- 
set at m„l2 provides an exact definition of the pion mass which 
is manifestly dififerent from the screening mass as pointed out in 
[44] already. The (square of the) screening mass is the inverse 
of the propagator at - 0. It diff'ers from the physical mass 
by the amount by which its radiative corrections change when 
extrapolating from p^ = to the pion pole at -p^ - m^. This 
is all very-well known, of course. What might be surprising, 
however, is that with the independent exact mass definition at 
hand, via the quantum phase transition, one observes that this 
difiference can be as much as 30% in the QM models [44]. The 
approximation to mistake the pion screening mass for the phys- 
ical one thus has a considerable Silver Blaze problem. 

In general, for more realistic model parameters, indepen- 
dent of the existence of a quantum phase transition, one there- 
fore needs to consider the poles in the zero-temperature prop- 
agators or, equivalently, the zeros in the corresponding 2-point 
functions, to define the physical masses. They are determined 
from the asymptotic behavior of the propagators at large Eu- 
clidean times. For the lowest stable particle it suffices to solve 
analytically continued flow equations at timelike total momenta, 
and find, e.g., for the pion. 



rf'^\p;cr^,,) 



= 0. 



(14) 



At mean-field level, the Silver Blaze problem requires this pion 
pole mass to be computed from the random phase approxima- 
tion (RPA). It is then guaranteed to agree with the critical chem- 
ical potential at the BEC quantum phase transition in NIL and 
QM models [44, 49]. Here, we use the ERG to calculate a pion 
pole mass in a way consistent with the derivative expansion of 
the effective average action. Although an analogous flow equa- 
tion can be obtained for the 2-point function of the sigma, a 
pole mass determination requires more work in that case and 
will not be done here. The additional complication is associ- 
ated with the finite width of the sigma meson from its decay 
into two pions, which requires a more careful search for the 
sigma pole on the correct sheet in the complex plane. 

Applying two functional derivatives to Eq. (6), one obtains 
the flow equations for the 2-point functions as sketched in Fig. 2, 
which depend on scale-dependent 3- and 4-point functions with 
their own flow equations. Truncations are necessary in order to 
close this infinite tower of equations for all «-point functions. 
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Figure 2: Diagrammatic representation of the flow equation for a mesonic 2- 
point function. Solid (dashed) lines represent scale-dependent fermion (boson) 
propagators and the light circles the insertions of d,Rk{q)- 



One example of such a truncation scheme is the so-called 
BMW approximation [50, 51]. The basic idea is that, because 
of the insertion of the regulator function, the dependence of 
scale-dependent 3- and 4-point functions on the external mo- 
menta is weaker than on the loop momentum. This motivates 
to expand 3- and 4-point functions in their external momenta. 
At leading order this yields, e.g., for the mesonic 3- and 4-point 
functions, 

(0 3) dTf^iP^d^ (0 4) (^^\P^1» 



d(p,„d(pi 



The resulting closed coupled system of flow equations for the 
2-point correlation functions has been studied for example for 
scalar models in [52, 53]. Consistency between the BMW and 
LPA approximations for vanishing external momentum, where 
the 2-point function can also be obtained from derivatives of the 
effective potential, is not guaranteed, however 

Therefore, we use an even simpler truncation which is man- 
ifestly consistent, however, with the leading order in the deriva- 
tive expansion for the effective average action. This is achieved 
by using the scale dependent but momentum independent me- 
sonic 3- and 4-point vertices from the computation of the scale 
dependent effective potential in the flow for the 2-point func- 
tions, i.e.. 
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With constant Yukawa coupling at the leading order deriva- 
tive expansion the quark-meson vertices take their simple RPA 
forms as in the NJL model [54], 



r(2,i) 
^ 



-(2,1) 



p(2,2) ^ Q 



(16) 



The momentum dependence of the 2-point correlation function 
now arises entirely from the external momentum through the 
propagators in the loop. For the vacuum pion mass at zero 
isospin and baryon chemical potential the effective potential 
remains (9(4) symmetric. We therefore replace Uk{p^,cf') by 
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Uki(/>^) with - + cfi in our Ansatz for the effective aver- 
age action. The LPA-Uke flow of the pion 2-point function for 
Euclidean external momentum p = (po, 0) then becomes. 
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with standard primes for the 0^-derivatives here in the place 
of the and t/^ -derivatives of Uk(p^,cfi) above, and E-„ = 
y/k^ + 2U', Ea- = yjk^ + 2U' +4(f>^U" correspondingly. We 
use N - A for the bosonic (9(4) symmetry with Nj - 2 flavors 
(as compared to N - 6 with Nf = 2 in two-color QCD [44]). 
At zero pion momentum one easily verifies that Eq. (17) agrees 
with the flow equation for 2f/^ obtained from Eq. (11). 

We solve the flow equation for the real-time 2-point func- 
tion by analytically continuing -/jy - Luf^ + is in Eq. (17) before 
the integration of the flow equation. The initial condition in the 
UV is then chosen to be 



= -o/ + 2U\ 



(18) 



3. 3. Extended mean-field ( eMF) approximation 



A worthwhile further simplification of the flow equations 
for the effective potential in Eq. (11) and the 2-point function 
in Eq. (17) is to consider the purely fermionic flow, i.e., to ne- 
glect the mesonic contributions to the flow. This is occasionally 
being referred to as extended mean-field approximation (eMF) 
in the ERG literature because the fermionic flow is independent 
of the mesonic potential. Moreover, it provides an unambigu- 
ous way to include vacuum contributions which affect the usual 
mean-field studies [44, 55]. The corresponding flow equations 
can then formally be integrated, yielding for the effective poten- 
tial at finite temperature, isospin and quark chemical potential. 
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Figure 3: Vacuum pion 2-point correlation function from tlie full FRG Eq. (17) 
as a function of timelike a) = ipo compared to i,f^i,a (~''^) ■ 



the screening mass mf 



V2l7', both from the FRG parameter set C in Tab. 3. 



4. Results 

We solve the flow equations Eq. (11) and Eq. (17) numer- 
ically on grids in field space. For example, in the (9(4) sym- 
metric case for Uk{<p^) and rf'^\p; (p) this is a one-dimensional 
set of discrete field values of <p^. The necessary derivatives of 
the effective potential can be obtained from finite differences or 
spline interpolations. This leads to a system of ordinary differ- 
ential equations which can be integrated with standard methods. 
In comparison to using Taylor expansions of the effective poten- 
tial around a scale dependent minimum, it captures the form of 
the effective potential everywhere on the grid in field space in- 
stead of being restricted to some region around the minimum. It 
is thus well suited to describe first order phase transitions. Only 
recently this technique has been applied to problems with effec- 
tive potentials depending on several invariants [44] . We use the 
same techniques here when we study the full two-dimensional 
flow from Eq. (11) for Ukip^, d^) with collective fluctuations in 
both, the chiral and the charged pion condensate. 

4.1. Pole mass versus screening mass 

Fig. 3 shows the square root of the modulus of the pion cor- 
relation function evaluated at the minimum of the effective po- 
tential cr,„in as a function of timelike external u - ip^. For com- 
parison, the corresponding T^^l^ = -of' + with the screen- 
ing mass m^'^'' from the curvature of effective potential is also 



ftanh ^ ^ ^ + tanh shown. Pole and screening mass are the zeros of rJr'^(-/(w; cr^in) 
V \ 2T ) \ 2T )/' and rger(-iw;o-„in). Moreover, from Eqs. (11) and (17), TO^"2 = 



Ua = a(p^ +d^) + b [p^ + d^) -co-- 2^i]d^ 



rScr(0; o-mm) = '^'^n'i^', cTmin), in our truncation. 

In the present example, the pole mass is m^^^^ - 133.0 MeV 
which is within 3% of the correct value as defined by the quan- 
tum phase transition at 2yuJ = 136.6 MeV with the parameters 
used here. The screening mass, on the other hand, is much 
heavier, with nr"" - 188.0 MeV it overestimates the correct 
value by about 38%. The Silver Blaze relation m„ - 2p'j is sat- 
isfied exactly for the pion pole mass in the eMF calculation, as it 
must [44]. For a full FRG calculation, the pion pole mass in our 
truncation is close and represents a considerable improvement 
as compared to the screening mass. Perfect agreement should 
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and for the vacuum pion 2-point function 

,2 



Ff 2>(-/w;(r) = dkk 



4El + 



E,{^El- 



-(2) 

' A,7T 



-a? + 2a + Abcr^. 



- + r^'' 

i2 A,;r' 



(20) 



As before, the real-time result is obtained via analytic continua- 
tion with -py = cxP'+iE from the Euclidean correlation function. 



not be expected in this calculation for the pion pole mass either, 
since the truncation for the flow of the 2-point function goes be- 
yond the leading order derivative expansion. Stricktly speaking 
one would have to feed the momentum dependent propagators 
back into the flow equation for the effective potential and obtain 
both in an iterative procedure. Maybe more importantly, how- 
ever, these results once more demonstrate that assigning phys- 
ical values to screening masses is a rather poor way of fixing 
model parameters. This is even more so for the pion mass with 
3 colors here than in the corresponding 2-color calculation [44]. 

Less pronounced, the same inconsistency between screen- 
ing mass and quantum phase transition at fi - ffj has been ob- 
served already in a purely bosonic model with isospin chemi- 
cal potential [56]. For a screening mass of m^" = 138 MeV, 
the zero-temperature onset of pion condensation is observed in 
this model at = 132.6 MeV, corresponding to a moder- 
ate mismatch of 4%. For comparison, we obtain a pion pole 
mass of mjj"''^ - 129.9 MeV in our truncation for this purely 
bosonic model. The resulting deviation of only 2% again rep- 
resents at least some improvement as compared to the screen- 
ing mass. The relatively small deviation of pole and screening 
masses here is in line with the general observation that fermions 
provide the dominant contributions to the flow in QM models. 

4.2. Parameter fixing 

eMF calculation. For the extended mean-field calculation, we 
select parameter sets A, B, C and D which reproduce the pion 
decay constant f„ and the pion pole mass and in the vacuum 
from Eqs. (19) and (20), as summarized in Tab. 2. In this range 
of parameters the mass m^r of the sigma meson increases mono- 
tonically with the quartic coupling b in the UV potential. 





fl/A^ 


b 


cjh? 


fn 


pole 

mn 




A 


0.245 


1.0 


0.0124 


92.8 


138 


457 


B 


0.196 


2.0 


0.0123 


92.9 


137 


504 


C 


0.09 


4.0 


0.0125 


92.4 


138 


698 


D 


-0.09 


8.0 


0.0125 


92.8 


139 


868 



Table 2: Parameters for eMF calculations with UV cutoff A = 600 MeV and 
Yukawa coupling g = 3.2; and meson masses are given in MeV. 

The left panel in Fig. 4 shows the resulting (jj., T) phase di- 
agram with Hi - and no pion condensation. For the param- 
eter sets A and B there is a critical endpoint (CEP) at (r,ju) = 
(22.5 MeV, 314 MeV) and (32 MeV, 295 MeV), respectively. 
The position of the CEP is very sensitive to the choice of pa- 
rameters. With increasing nicr the CEP moves to lower temper- 
atures and eventually disappears from the phase diagram be- 
yond a critical value of m^"* ~ 550 MeV, so there is none for 
parameter sets C and D. A similar behavior was observed in a 
mean-field 3-flavor QM model calculation [57]. An intuitive 
understanding of this effect is gained by considering a Landau 
expansion of the effective potential in cr^, yielding, 

U((r, T,i^) = a'(T,H)(r^ + b'iT,^i)cr' + d(T,i^)a-\ 
a'(T, ij)^a(T,iJ.) + a, b'(T, /j.) ^b{T,n) + b. 




H[Ml-V] ^i(MeV| 



Figure 4: QM model phase diagram at ,u; = from eMF (left, parameter sets 
A,B,C, D from left to right) and full FRG (right, parameter sets A,B,C) calcu- 
lations. Solid lines represent first order phase boundaries and dashed lines trace 
minima in the screening mass of the sigma meson to indicate chiral crossovers. 

To describe a first order phase transition we need to keep terms 
up to sixth order The corresponding quark contributions are 
a{T,ii), b(T,fi) and d{T,fj.). They are independent of the initial 
parameter choices for a and b. Around the first order phase 
boundary, the total coefficients satisfy a' > 0, b' < and d > 0, 
which is just the condition for a double-well. For a sufficiently 
large quartic bare coupling b, this condition can no longer be 
satisfied. Thus the first order phase transition and the CEP have 
to disappear from the (/u, T) phase diagram as m,j- increases. 

FRG calculation. Three sets for initial parameters in full FRG 
solutions are listed in Tab. 3. As for the eMF results, the screen- 
ing mass of the sigma monotonically increases with the quartic 
coupling parameter b. 





a/A^ 


b 


c/A^ 




pole 

K 


K' 


A 


0.3224 


0.25 


0.0045 


93.3 


134.3 


524.1 


B 


0.2844 


1.25 


0.0045 


92.8 


134.1 


557.1 


C 





9.075 


0.0045 


92.8 


133 


696.7 


A 


-0.48 


9.74 


0.0024 


93.0 


138.0 


448 



Table 3: Initial parameters for full FRG calculations (parameter set A' for the 
purely bosonic model). The UV cutoff is set to A = 900 MeV and g = 3.2; 
and meson masses are given in MeV. 



The right panel in Fig. 4 shows the yU/ = phase diagram 
in the (/i, T) plane in the vicinity of the CEP for the three pa- 
rameter sets A,B and C. In contrast to the eMF calculation, the 
phase structure for all sets is quite similar, although the screen- 
ing masses are rather different. In all cases there are first order 
phase transitions and CEPs. In the full FRG effective potential, 
it makes no sense to separate the initial potential from the con- 
tributions due to fluctuations as in Eq. (21). Differences at the 
UV scale are reduced by the mesonic fluctuations. As a result, 
the phase structure from full FRG calculations is less sensitive 
to the resulting sigma mass than in mean-field calculations. 

4.3. Zero temperature phase diagram 

First consider the zero temperature quantum phase transi- 
tion with pion condensation ai ji - 0. We have calculated the 
isospin density from the ///-derivative of the effective potential. 
The result is shown in Fig. 5 where the rescaled isospin chem- 
ical potential 2iij/m„ is plotted over the isospin density pj for 
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Figure 5: Isospin chemical potential over isospin density at zero temperature 
and baryon cliemical potential in comparison to the lattice data from Ref. [17]. 
The isospin densities are rescaled here to adjust to the lattice parameters. This 
mainly compensates for the larger pion mass and decay constant on the lattice. 

our eMF and FRG solutions in comparison with very recent lat- 
tice data from Ref. [17] and the leading order chiral effective 
field theory O^fPT) result [12]. The isopin density remains zero 
below the onset of pion condensation at jj'j - m„l2 in the eMF 
calculation as it must while the full FRG result shows a slight 
residual Silver Blaze problem here, which might be due to nu- 
merical integration and infrared cutoff uncertainties. For large 
yU/, deep into the pion condensation phase, it gets increasingly 
difficult to control systematic errors in the full FRG calcula- 
tions due to, e.g., an increasing sensitivity to the sigma mass 
assignment and the ultraviolet cutoff. 

While the chiral effective field theory result describes the 
BEC phase well, it appears to miss some essential dynamics 
in the BEC-BCS crossover, which we estimate to occur around 
2fj.j/m„ - 1 X! 2/3 in our calculations by the simple criterion 
that the quarks' Dirac mass falls below /i/ there. Because the 
vacuum realignment in this crossover is characterized by the 
mixing between the sigma meson with the charged pions in the 
QM model, this suggests that it might indeed be this mixing 
of the meson mass eigenstates which is responsible for the ob- 
served inflection point and back-bending of the isospin density 
at large fij. The analogous mixing is seen in the mass spectrum 
of the QMD model for two-color QCD [44], between diquarks 
and the sigma meson. Basically, the latter is infinitely heavy 
in ;i'PT, so the low-lying excitation spectrum of the QM model 
is essentially different from that of the non-linear sigma model, 
especially in the crossover region where the mixing occurs. The 
qualitative effect on the isospin density can be illustrated in the 
linear sigma model, where a simple calculation yields, 

p/(jc, y) = 2f„ m„ X [yZi 'p^^ I ' ^^^"^ 

for X = Ijiilnij, > 1, and with y = ma-/m„. In the y — > oo limit, 
this reduces to the;i'PT result pi(x) - 2f^m„x{\ - jc"'*) [12]. 
Compared to that, a finite sigma mass y < oo in the linear sigma 
model has the same qualitative effect as observed in the eMF 
results of Fig. 5 for parameter sets A,B and C in Tab. 2. 



The general structure of the phase diagram in the T - Q 
plane is also strongly constrained by the Silver Blaze property. 
The key to its understanding is the fact that different degrees of 
freedom couple to different combinations of the chemical po- 
tentials. While up(down)-quarks couple to /i + /i/, charged pi- 
ons obviously only couple to the isospin chemical potential /i/. 
The partition function and, correspondingly, thermodynamic 
observables must remain independent of and fij as long as 
jj + fdi < mg = go-^ia = gf„ and fij < m„/2, where nig and 
m„ are the vacuum quark and pion masses, respectively. This 
defines a quadrilateral area in the zero-temperature phase dia- 
gram, which we refer to as first Silver Blaze region in the fol- 
lowing. The physical picture is that in absence of bound quark 
matter, the horizontal line in the (jU/, //) plane defines the bound- 
ary of this region beyond which a degenerate Fermi gas of up- 
quarks forms, just as the vertical line marks the onset of pion 
condensation. Another constraint arises inside the pion conden- 
sation phase: For /// > m,r/2 the vacuum changes and the grand 
potential becomes ///-dependent, of course. Nevertheless, for 
constant jjj, it has to remain independent of the isosymmetric 
quark chemical potential fi as long as this stays below the light- 
est quark mass, i.e., for fi < nigi/xi), where 

m^gifi,) = ylg^d^ + (gp±fi,)2 (23) 

are the quark masses at yu = 0. They no longer correspond 
to pure up- and down- quark excitations as these mix in the 
pion condensation phase, and the quark masses are obtained by 
diagonalizing r^^'*'\ More generally, from the quark dispersion 
relation 



Egif) = ^g'd^ + ( ^p'+gV - /^/j , (24) 

one furthermore notes that for /// > gp it becomes energeti- 
cally favorable to excite quarks with a finite spatial momentum 

= yu^ - g^p^ and minimal energy ■ - gd [37]. This 
might be an indication for an inhomogeneous phase at large 
isospin chemical potential such as the FFLO phases [27, 28] 
also discussed for two-color QCD with isospin chemical poten- 
tial [58, 59] or in the 1 + 1 dimensional NJL model [60]. 

Therefore, the conditions jii > ni„/2 with fi < nig(fij) for 
fii < gp, or yu < E'. for > gp, together define a region in 
the zero-temperature phase diagram in which the ///-dependent 
grand potential still remains independent of //, nevertheless. We 
will call this the second Silver Blaze region in the following. 

Note however that the whole argument holds only as long 
as there are no first order transitions in these regions of the 
phase diagram. This concerns in particular the possible first 
order transition for large fi and small /// which may intersect 
with the first Silver Blaze region. If it does, the CEP where it 
ends must lie outside, however This could then be the chiral 
first order transition or a liquid-gas transition to bound quark 
matter below the threshold for free quarks, for example, which 
are of course not excluded. In contrast, a first order line with 
an endpoint inside the Silver Blaze region would be thermody- 
namically inconsistent. It is reassuring that we never observe 
this in our numerical calculations either 
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Figure 6: (fii,/!) phase diagram at T = from an eMF calculation (parameter 
set A). Tlie SB line marks the boundary of the first/second Silver Blaze region. 
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Figure 7: (jJi,fi) phase diagram at T = from the full FRG calculation (param- 
eter set C). The SB line represents fi = m^. 



eMF calculation. Fig. 6 shows the phase diagram in the (fi[,jJ.) 
plane from an eMF calculation. As discussed above, there is a 
chiral first order transition outside the pion condensation phase 
for large baryon chemical potential and small isospin chemi- 
cal potential. With parameter set A, the first order transition 
Ues completely outside the first Silver Blaze region bounded by 
1-1+ fii < nig. The vertical line in the phase diagram separates the 
charged pion condensation phase from the normal chiral sym- 
metry breaking phase. This phase boundary stays at the con- 
stant fij = ni„/2 until fi > nig - m„/2 from where on it bends 
to larger values of /i/ when further increasing jj. For fi; > in„/2 
the general arguments from above require the phase boundary 
of the pion condensation phase, as long as it is of second or- 
der, to stay outside the second Silver Blaze region. Indeed, in 
perfect agreement with these general arguments, the boundary 
of the pion condensation phase crosses the Silver Blaze line 
only after it has become of first order, beyond a tricritical point 
(TCP) which is here observed at = (0.56 m;,, 253 Me V). 

We will see, however, that this first order line is washed out by 
mesonic fluctuations, and that it thus is a mean-field artifact. 

At the mean-field level the Silver-Blaze property is manifest 
in the explicit expressions for the partition function. This has 
been observed in [41] as a property of the zero temperature par- 
tition function. It also holds for the explicit eMF expression for 
the effective potential Eq. (19): one verifies that it remains con- 
stant throughout the first Silver Blaze region with /i + /i/ < nig 
and fij < m„. Moreover, with Eq. (19), the condition for the 
second-order phase boundary of charged pion condensation. 



0, 



(25) 



which defines the critical isospin chemical potential jfjiii), is in 
fact independent of fi for/i + fii < gf„ = nig. This implies that 
the boundary of the pion condensation phase stays constant at 
fJ^i = iii„/2 for yu < nig - ni„/2 as also seen in Fig. 6. In fact, 
Eq. (25) coincides with 



(26) 



from Eq. (20) there. This establishes explicitly that the pole 
mass agrees with the onset of pion condensation at f/j = ni„/2 
for all ju < nig - ni„l2 in the eMF calculation. 

FRG calculation. Fig. 7 shows the corresponding phase dia- 
gram from the full FRG calculation. As in mean-field calcu- 
lations there is a first order transition at large ^ and small ///. 
However, the mesonic fluctuations tend to weaken this first or- 
der transition and correspondingly the CEP is at a smaller 
than in eMF. For small chemical potential ji, the second order 
transition separating charged pion condensation from the nor- 
mal ;^'SB phase again occurs at yt// = ni„l2. And again, the 
phase boundary stays constant until it hits the first Silver Blaze 
boundary line at ju = nig-ni„/2, which defines the edge between 
the first and the second Silver Blaze region. 

In contrast to mean-field calculations, however, the pion 
condensation transition with the mesonic fluctuations remains 
second order for the whole range of isospin chemical potentials 
we have investigated. Note also that this second order phase 
transition occurs outside the Silver Blaze region. Interestingly, 
one observes another first order transition inside this region for 
large isospin chemical potentials which was not present in the 
eMF calculation. This implies that the zero temperature parti- 
tion function and correspondingly also the chiral and charged 
pion condensates at fixed /// stay independent of /i only below 
this first order transition line which seems to end in a CEP right 
at the boundary of the second Silver Blaze region. 

4.4. Three dimensional T) phase diagram 

At finite temperature the phase boundaries in Figs. (6) and 
(7) extend into surfaces in the three dimensional parameter space 
with temperature T over the (ji!,fi) plane. These are shown for 
a representative eMF calculation in Fig. 8, and a full FRG solu- 
tion on the two-dimenional grid in field space with fluctuations 
in the chiral as well as the charged pion condensate in Fig. 9. 

The green surfaces in these figures represent the boundary 
of the pion condensation phase. The zero temperature TCP of 
Fig. 6 in the eMF calculation extends into the purple line in 
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Figure 8: Three dimensional pliase diagram from eMF (parameter set A). 



Figure 9: Three dimensional phase diagram from full FRG (parameter set C). 



Fig. 8, which divides the surface of the pion condensation tran- 
sition into a second order region (dark green, smaller yU/) and a 
first order region (light green, larger At the other end of this 
purple tricitical line, in the ipi, T) plane at = 0, it becomes the 
analogue of a TCP that was also predicted for two-color QCD 
at finite baryon chemical potential from next-to-leading order 
XPT [61]. As in the QMD model for two-color QCD, where 
the fluctuations of collective mesonic and baryonic excitations 
wash out the corresponding first order line [44], we observe 
that the chiral and charged pion condensate fluctuations have 
the same effect here. In the QM model for QCD with baryon 
and isospin chemical potentials this happens for all values of yu, 
starting from the T = end to the two-color like yu = plane. 
The whole (light green) first order surface towards the larger jjj 
values is gone in Fig. 9 with fluctuations, and with it the entire 
tricritical line. 

The chiral first order line of the T - eMF phase diagram 
for small ///, just outside the first Silver Blaze region in Fig. 6, 
is the baseline of the red surface in the corresponding three- 
dimensional eMF phase diagram in Fig. 8. Its edge represents 
the line of the critical endpoint in Fig. 4 which moves to lower 
temperature with increasing fij ending up as that in Fig. 6. As 
one can also see from these figures, it thereby stays well clear 
of the pion condensation phase. The corresponding (red) first 
order surface with fluctuations in Fig. 9 by and large shows a 
similar behavior, albeit being much smaller than in the eMF 
calculation as already observed atT -0. This first order transi- 
tion does not disappear but it is considerably weakened by the 
fluctuations as well. 

With fluctuations, there is another first order surface inside 
the pion condensation phase. This is the extension into the three 
dimensional parameter space with temperature of the first or- 
der line near the boundary of the second Silver Blaze region in 
Fig. 7. We have not investigated it in much detail here and it is 
not included in Fig. 9, but it rises only to relatively small tem- 
peratures ending in another line of critical endpoints similar to 
that of chiral first order transition, but entirely contained inside 
the pion condensation phase. 

Finally, in Fig. 10 we show various (yU/, T) planes for differ- 



ent values of the isosymmetric quark chemical potential yU of the 
FRG phase diagram with fluctuations in Fig. 9. The fi -Q plane 
thereby corresponds to the (jj, T) phase diagram of the QMD 
model for two-color QCD with diquark instead of charged pion 
condensation. This naturally reflects the precise map between 
the corresponding flow equations. Fig. 10 also shows how the 
chiral condensate vanishes as the vacuum rotates in the pion 
condensation phase. For small baryon chemical potentials, be- 
low the edge of the Silver Blaze regions at yU = mq-m„/2, small 
quantitative changes arise at finite T. At larger the imbalance 
between the Fermi surfaces of up and anti-down quarks builds 
up, and the pion condensation phase starts to shrink rapidly, 
likely leading into an FFLO phase or alike in some possibly 
smaU region where both fi and fi; are about equally large. 
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Figure 10: (p/,T) slices for different baryon chemical potentials of the three- 
dimensional FRG phase diagram in Fig. 9. Also shown here are chiral crossover 
lines (dashed) as the half-value of the chiral condensate. 



5. Summary and Conclusions 

We have investigated the phase diagram of the two-flavor 
QM model for QCD with both isospin and baryon chemical 
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potential using the FRG. We showed explicitly how this effec- 
tive model can be mapped onto the QMD model for two-color 
QCD. As compared to earlier NJL model studies, mesonic fluc- 
tuations are included with the FRG. We have demonstrated their 
efl'ects on the phase diagram by detailed comparisons of purely 
fermionic flows in the eMF approximation with full FRG re- 
sults. As compared to previous FRG studies, collective meso- 
nic excitations are included in a way suitable to simultaneously 
describe the competing fluctuations in both, the chiral and the 
charged pion condensates, based on recent developments within 
the QMD model for two-color QCD. As compared to the latter, 
a finite baryon chemical potential here adds another dimension 
to the parameter space of the phase diagram. Without baryon 
chemical potential, our results compare well with recent lattice 
data indicating that the sigma might be important for the me- 
son mixing and the change in the low-lying excitation spectrum 
over the BEC-BCS crossover in the pion condensation phase. 
At zero temperature, baryon chemical potentials below gener- 
alized Silver Blaze bounds have no efl'ect. Beyond those, the 
increasing imbalance of the up and anti-down Fermi energies 
reduces pion condensation as expected. Moreover, the onset of 
pion condensation once more demonstrates how important it is 
to use proper pole masses rather than screening masses for a 
realistic parameter fixing in QM models. We have shown how 
to obtain good estimates for those from the FRG in a way that 
is consistent with the FRG calculation of the grand potential. 

Acknowledgements 

K.K. was supported by the Grant-in-Aid for ISPS Fellows (No. 
22-3671) and by fliat for the Global COE Program "The Next 
Generation of Physics, Spun from Universality and Emergence." 
This work was supported by the Helmholtz International Cen- 
ter for FAIR within the LOEWE program of the State of Hesse, 
the Helmholtz Association Grant VH-NG-332, and the Euro- 
pean Commission, FP7-PEOPLE-2009-RG No. 249203. 

References 

[1] p. Braun-Munzinger and J. Wambach, Rev. Mod. Phys. 81 (2009) 1031. 
[2] B. Friman, C. Hohne, J. Knoll, S. Leupold, J. Randrup, R. Rapp and 

P Senger, Lect. Notes Phys. 814 (2011) 1. 
[3] K. Fukushima and T. Hatsuda, Rept. Prog. Phys. 74 (2011) 014001. 
[4] Ph.deForcrand,PoSLAT2009 (2009)010. 
[5] M. Buballa, Phys. Rept. 407 (2005) 205. 
[6] K. Fukushima, Phys. Lett. B 591 (2004) 277. 

[7] C. Ratti, M. A. Thaler and W. Weise, Phys. Rev. D 73 (2006) 014019. 

[8] J. Herges, D. Jungnickel and C. Wetterich, Eur. Phys. J. C 13 (2000) 323. 

[9] J. Braun, H.-J. Pimer and K. Schwenzer, Phys. Rev. D 70 (2004) 085016. 
[10] B.-J. Schaefer and J. Wambach, Nucl. Phys. A 757 (2005) 479. 
[11] B.-J. Schaefer, J. M. Pawlowski and J. Wambach, Phys. Rev. D 76 (2007) 
074023. 

[12] D. T. Son and M. A. Stephanov, Phys. Rev. Lett. 86, 592 (2001). 

[13] J. B. Kogut, M. A. Stephanov, D. Toublan, J. J. M. Verbaarschot and 

A. Zhitnitsky, Nucl. Phys. B 582 (2000) 477. 
[14] J. B. Kogut and D. K. Sinclair, Phys. Rev. D 70 (2004) 094501. 
[15] D. K. Sinclair and J. B. Kogut, PoS LAT 2006 (2006) 147. 
[16] Ph. de Forcrand, M. A. Stephanov and U. Wenger, PoS LAT 2007 (2007) 

237. 

[17] W. Dermoid, K. Orginos and Z. Shi, arXiv: 1205.4224 [hep-latj. 
[18] A. Nakamura, Phys. Lett. B 149 (1984) 391. 



[19] S. Hands, J. B. Kogut, M. P. Lombardo and S. E. Morrison, Nucl. Phys. 
B 558 (1999) 327. 

[20] S. Hands, et at, Eur. Phys. J. C17 (2000) 285-302; S. Hands, I. Montvay, 
L. Scorzato, and J. Skullerud, Eur. Phys. J. C22 (2001) 451-461. 

[21] S. Hands, R Kenny, S. Kim and J. I. Skullerud, Eur. Phys. J. A 47 (201 1) 
60. 

[22] K. Holland, R Minkowski, M. Pepe and U. J. Wiese, Nucl. Phys. B 668 

(2003) 207. 

[23] A. Maas, L. von Smekal, B. Wellegehausen and A. Wipf, 

arXiv: 1203.5653 [hep-lat]. 
[24] L. von Smekal, arXiv:1205.4205 [hep-ph]. 

[25] J. Braun, L. M. Haas, F. Marhauser and J. M. Pawlowski, Phys. Rev. Lett. 

106(2011)022002. 
[26] J. M. Pawlowski, AlP Conf. Proc. 1343 (201 1) 75. 
[27] M. G. Alford, J. A. Bowers and K. Rajagopal, Phys. Rev. D 63 (2001) 

074016. 

[28] P Fulde and R. A. Ferrell, Phys. Rev. 135 (1964) A550; A. I. Larkin and 
Y. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 47 (1964) 1 136 [Sov. Phys. JETP 
20(1965) 762]. 

[29] M. Hanada and N. Yamamoto, arXiv: 1103.5480 [hep-ph]. 

[30] D. T. Son and M. A. Stephanov, Phys. Atom. Nucl. 64 (2001) 834 [Yad. 

Fiz. 64 (2001) 899]. 
[31] M. C. Birse, T. D. Cohen and J. A. McGovem, Phys. Lett. B 516 (2001) 

27. 

[32] B. Klein, D. Toublan and J. J. M. Verbaarschot, Phys. Rev. D 68 (2003) 
014009. 

[33] T. Kanazawa, T. Wettig and N. Yamamoto, JHEP 1112 (201 1) 007. 
[34] T. D .Cohen, Phys. Rev. Lett. 91 (2003) 222001; arXiv;hep-ph/0405043. 
[35] A. Barducci, R. Casalbuoni, G. Pettini and L. RavagU, Phys. Rev. D 69 

(2004) 096004. 

[36] L. -y. He, M. Jin and P. -f. Zhuang, Phys. Rev. D 71 (2005) 116001. 
[37] D. Ebert and K. G. Klimenko, J. Phys. G G 32 (2006) 599. 
[38] Z. Zhang and Y. X. Liu, Phys. Rev. C 75 (2007) 064910. 
[39] S. Mukheijee, M. G. Mustafa and R. Ray, Phys. Rev. D 75 (2007) 094015. 
[40] J. Xiong, M. Jin and J. Li, J. Phys. G 36 (2009) 125005. 
[41] T. Sasaki, Y. Sakai, H. Kouno and M. Yahiro, Phys. Rev. D 82 (2010) 
116004. 

[42] C. -f. Mu, L. -y. He and Y. -x. Liu, Phys. Rev. D 82 (2010) 056006. 
[43] J. I. Kapusta and C. Gale, Cambridge, UK: Univ. Pr. (2006) 428 p. 
[44] N. Strodthoff, B.-J. Schaefer and L. von Smekal, Phys. Rev. D 85 (2012) 
074007. 

[45] C. Wetterich, Phys. Lett. B 301 (1993) 90. 

[46] J. Berges, N. Tetradis and C. Wetterich, Phys. Rept. 363 (2002) 223. 
[47] J. M. Pawlowski, Annals Phys. 322 (2007) 2831. 
[48] D. R Litim, Phys. Rev. D64 (2001) 105007. 

[49] T. Brauner, K. Fukushima and Y. Hidaka, Phys. Rev. D 80 (2009) 074035 

[Erratum-ibid. D 81 (2010) 119904]. 
[50] J. P Blaizot, R. Mendez-Galain and N. Wschebor, Phys. Rev. E 74 (2006) 

051116. 

[51] J. P. Blaizot, R. Mendez-Galain and N. Wschebor, 

Phys. Rev. E 74 (2006) 051117. 
[52] D. Guerra, R. Mendez-Galain and N. Wschebor, Eur. Phys. J. B 59 (2007) 

357. 

[53] F. Benitez, R. M. Galain and N. Wschebor, Phys. Rev. B 77 (2008) 
024431. 

[54] T. Hatsuda and T. Kunihiro, Phys. Rev. Lett. 55 (1985) 158. 

[55] V. Skokov, B. Friman, E. Nakano, K. Redlich and B. -J. Schaefer, Phys. 

Rev. D 82 (2010) 034029. 
[56] E. E. Svanes, J. O. Andersen, Nucl. Phys. A857 (201 1) 16. 
[57] B.-J. Schaefer and M. Wagner, Phys. Rev. D 79 (2009) 014018. 
[58] K. Fukushima and K. lida, Phys. Rev. D 76 (2007) 054004. 
[59] J. O. Andersen and T. Brauner, Phys. Rev. D 81 (2010) 096004. 
[60] D. Ebert, N. V. Gubina, K. G. Klimenko, S. G. Kurbanov and 

V. C. Zhukovsky, Phys. Rev. D 84 (2011) 025004. 
[61] K. Splittorff, D. Toublan and J. J. M. Verbaarschot, Nucl. Phys. B 639 

(2002) 524. 



10 



